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We study theoretically the entanglement created in a scattering between an electron, incoming 
from a source lead, and another electron bound in the ground state of a quantum dot, connected to 
two leads. We analyze the role played by the different kinds of resonances in the transmission spectra 
and by the number of scattering channels, into the amount of quantum correlations between the 
two identical carriers. It is shown that the entanglement between their energy states is not sensitive 
to the presence of Breit-Wigner resonances, while it presents a peculiar behavior in correspondence 
of Fano peaks: two close maxima separated by a minimum for a two-channel scattering, a single 
maximum for a multi-channel scattering. Such a behavior is ascribed to the different mechanisms 
characterizing the two types of resonances. Our results suggest that the production and detection of 
entanglement in quantum dot structures may be controlled by the manipulation of Fano resonances 
through external fields. 
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I. INTRODUCTION 

Quantum entanglement, as one of the most spectac- 
ular features of quantum mechanics contrasting with 
classical physics 1 -, has been widely investigated in the 
last decades, mainly because it is recognized as a cru- 
cial resource for quantum information processing and 
quantum communication^. It is therefore a problem of 
great interest to find physical systems where the entan- 
glement can be produced, manipulated and detected. 
Recently there have been different proposals to pro- 
duce entangled states, such as those based on atomic 
systems 3 , quantum electrodynamics cavities 4 - and solid 
state devices^iSiSiSiifl. Among the systems proposed for 
the realization of quantum information processing de- 
vices, those constructed with semiconductor quantum 
dots (QDs) are extremely promising due mainly to the 
controllability of their quantum stat o 11 ' 12 ' 13 . Indeed, 
semiconductor QDs posses many desirable features: an 
atomic-like structure, that can be fully controlled by an 
external electrostatical potential; a tunable coupling to 
source and drain leads, which makes feasible the integra- 
tion with other microelectronics devices; and the scala- 
bility, which seems to promise sophisticate engineering of 
the multi-QD structures. 

QDs represent also an ideal laboratory to compare 
exact numerical simulations of the quantum transport 
phenomena with experiments, since the number of de- 
grees of freedom involved is often small and the dis- 
creteness of QD states highly reduces the computational 
burden needed. Various effects such as the conductance 
quantization^ 4 -, the Coulomb blockade due to the elec- 
tron repulsion 1 ^, the interplay between resonances and 
the charging in QD structure d 16 ' 17 , strongly affect the 
transport properties. Another peculiar feature of the 
electron transmission through QDs is the partial reten- 



tion of quantum coherence 1 ^, whose measurement, by ex- 
perimental setups exploiting quantum interference (e.g. 
with a QD embedded in an Aharonov-Bohm interferome- 
ter), may yield information about transport phenomena, 
not readily available from conductance measure d 19 ' 20 ' 21 . 

In the frame of quantum transport in semiconductors, 
theoretical and experimental investigations have revealed 
other mechanisms governing the electronic transmission 
through QD structures. In particular it has been ob- 
served that two kinds of resonances can be present in 
conductance spectra, known as Fano and Breit-Wigner 
resonance a 22 ' 23 ' 24 ' 25 ' 26 . The formers are present when 
two transmission channels, a resonant one and a non- 
resonant one, interfere^ 7 -. Moreover they exhibit typical 
asymmetrical line shapes, with the transmitted phase in- 
creasing by 7r on the resonance peak and then dropping 
abruptly. The latters show a symmetric line shape and 
may be considered a limiting case of Fano resonances, 
occurring, for example, when there is a single channel 28 . 

In this paper we address the problem of entanglement 
generation in a two-particle scattering in a QD structure. 
We analyze, in particular, the role played by some mech- 
anisms of charge transport in the appearance of quan- 
tum correlations. To this aim we consider a scattering 
event in a one-dimensional (ID) double-barrier resonant 
tunnelling device (that mimics the confining potential of 
the QD), with an electron incoming from one lead and 
another electron bound in the ground state of the QD. 
The two particles feel the confining potential inside the 
device and interact through the Coulomb repulsion. In- 
deed carrier-carrier entanglement has been recently in- 
vestigated in various QD structures^ ' 10 ' 29 ' 30 ' 31 ' 32 ' 33 where 
different scattering setups are considered for the gener- 
ation of two-electron entangled states. In the present 
work, we adopt a time- independent few-particle approach 
that, although computational demanding, can be solved 
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numerically to obtain the exact modulus and phase of 
the transmission coefficient (TC) of an electron crossing 
the charged QD. This gives us the possibility to quan- 
tify the quantum correlation between the energy states 
of the scattered electron and the bound one, and to ex- 
pose its connection with the resonances exhibited by the 
TCs of the various scattering channels. Unlike previ- 
ous work s 8 ! 33 , the few-particle approach we use in this 
paper allows us to study the relation between the differ- 
ent kinds of resonances, or the number of energy levels 
available, with the entanglement formation. Even if the 
dynamics of carriers has been considered as ID, we can 
assume that the results obtained describe a general be- 
havior concerning also quantum transport in 2D and 3D 
physical systems, being the appearance of quantum cor- 
relations closely related to the nature of transport res- 
onances, whose underlying mechanisms are independent 
from the dimensionality of the system. 

The paper is organized as follow. In Sec. |TT] we de- 
scribe the physical model and the numerical approach 
adopted to calculate the two-particle scattering state and 
the quantum correlations in terms of the von Neumann 
entropy. In Sec. IIIII we present the numerical results ob- 
tained for the entanglement in the case of two differ- 
ent kinds of processes, namely two- and multi-channel 
scattering^ 3 -. Finally in Sec. IIVI we comment the results, 
draw final remarks and point out issues that require fur- 
ther research. 



II. THE PHYSICAL MODEL AND THE 
NUMERICAL APPROACH 

Our aim is to evaluate the entanglement between 
the energy states of two electrons, interacting via the 
Coulomb potential, one bound in a QD and the other 
passing through it. Here we summarize the physical 
model adopted for the open QD and the numerical ap- 
proach used to evaluate the two-particle scattering state 
and the entanglement. 

We consider a quasi ID double-barrier resonant tun- 
nelling device as, for example, the ones formed by ma- 
terial modulation in a ultrathin cylindrical nanowirc 34 . 
The transversal dimensions of the structure are small 
compared to the other lenghtscales so that a single 
transversal subband is accessible to the carriers and the 
effective dynamics can be considered ID 35 . Two small 
potential barriers separate the QD region from the two 
contacts, as depicted in Fig. [lja). The N bound states 
and energies of the QD will be indicated as Xn and E n 
(with n = 0, 1, . . . , N in order of increasing energy), re- 
spectively. 

A single electron is in the QD ground state xo, whereas 
a second electron is incoming from the left lead with en- 
ergy Ej 7v , and it is scattered by the structure potential 
V s and by the Coulomb interaction with the bound parti- 
cle. The potential V s is supposed to be constant outside 
the region of interest of length A and, without loss of 



generality, it is taken to be zero in both left and right 
contacts. We will consider only cases in which the en- 
ergy of the incoming electron is not sufficient to ionize 
the QD, i.e. (Ejn + Eq) < 0. This means that when an 
electron leaves the scattering region, either reflected or 
transmitted, the other one is in a bound state of the dot. 
The two-particle Hamiltonian reads: 



H(xi,x 2 ) 



2m 

+V s ( Xl ) + V s (x 2 ) + 



2m* 



(1) 



47reA/ (xi - x 2 ) 



d 2 



where m* and e are the electron effective mass and the di- 
electric constant of the material, respectively. In particu- 
lar the calculations presented in this paper have been per- 
formed using GaAs material parameters. The Coulomb 
term includes the thickness a cut-off term d that can be 
assumed to correspond roughly to the lateral dimension 
of the confinement^. The spin degree of freedom does 
not enter into our calculation since we neglect spin-orbit 
coupling 44 . As a consequence of the fermionic nature 
of the system we impose antisymmetrization constraints 
under particle exchange to the wavefunction ip(xi,x 2 ). 
This is done by adopting antisymmetric boundary con- 
ditions as briefly described in the following and detailed 
in Ref. M- 

The two-particle scattering state is obtained by solving 
the time-independent open-boundary Schrodinger equa- 
tion Hi/j = Eip in the 2D domain of interest. To this aim 
we have used a numerical approach, based on a general- 
ization of the widely used quantum transmitting boundary 
method^. It allows us to include proper open boundary 
conditions and simulate the scattering of one electron by 
a charge confined in a QD 25 . In particular here we need 
four boundary conditions, one for each edge of the square 
domain. Since we have to impose exchange symmetry to 
ip(xi, x 2 ), they are equal in couples, apart from the sign. 
In fact the form of the wavefunction when particle 1 is 
in the left lead (x\ < 0) is 

(2) 



ip(x 1 ,x 2 )\ {xi< Q ) = Xo(x 2 y koXl + 



M 

£ 

71=0 



bnXn{x 2 )e 



-ik n x 



1 b nXn(x 2 )e k ^, 



and, when it is particle 2 the one in the left 
lead, the boundary condition is t/j(xi, x 2 )\^ X2<0 - ) — 
—ip(x 2 , £i)|(x 2 <0)- For the other two boundaries, namely 
the conditions (xi > A) and (x 2 > A) it holds 

^(X!,X 2 )\ {X1>A) = 



M 

E 

n=0 



CnXn(x 2 )e lk " Xl + C nX n(x2)e- k " Xl (3) 

n=M+l 



and i>(xi,x 2 )\( X2> A) = -tp(x2, xi)\(x 2 >A), respectively. 

Let us now describe the above expressions. We first 
define T n = E IN + E — E n as the energy of an elec- 
tron freely propagating in the lead when the other one 
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is bound in Xn (a consequence of energy conservation), 
and k„ = v / 2m*\T n \/h. The right hand side of Eq. © 
(for the left boundary) is the sum of three terms. The 
first one represents an electron incoming as a plane wave 
with energy Ein = ft 2 fcg/2r7i*, while the other electron 
is in the QD ground state xo- The second term accounts 
for all the energy-allowed possibilities with one electron 
bound in the Xn state, and the other one reflected in the 
left lead, with wave vector k n . M is the number of states 
for which T n is positive. The third term accounts for the 
cases T n < 0, representing the electron in the lead as 
an evanescent wave. The right hand side of Eq. ((3]) (for 
the right boundary) has only two terms, since the prob- 
ability amplitude of a carrier incoming from the right 
lead is zero in our system. The first term represents the 
M energy- allowed possibilities of an electron transmitted 
and freely propagating in the right lead and the second 
term includes the outgoing-particle evanescent waves, as 
in Eq. ([2]). The reflection and transmission amplitudes 
in the various energy levels, b n and c„ respectively, are 
unknowns and are obtained by solving the Schrodinger 
equation with H(xi,X2) given by Eq. (TTJ) and with the 
two-particle energy E = Ein + Eq, imposing the two 
boundary conditions of Eqs. and ©. 

In order to describe the mechanisms characterizing the 
charge transport through a QD, and the kind of reso- 
nances showing up in transmission spectra, we report in 
Fig. [T] the modulus and phase of the TC as a function 
of the initial energy of the incoming electron, for a given 
configuration of the potential sketched in the panel (a). 
In particular we consider a width L of the potential well 
of 80 nm with a depth h of 150 meV. In this case we have 
a single-channel scattering, i.e. the scattered particle can 
have a single energy, while the bound one is always left 
in the ground state of the QD well. When Ein is around 
1.5 meV, the TC shows a Breit-Wigner resonance, as 
shown in panel (b)2£. As it is well known, such kinds 
of resonances stem from the coupling of a quasi-bound 
state to the scattering states in the leads and presents a 
Lorentzian line-shape whose amplitude Co is described 
by the expression c = C{iT /2)/(E IN — E qb + iT/2), 
where C is a complex constant, T the width of the reso- 
nance, inversely proportional to the lifetime of the quasi- 
bound state with energy E q b. In particular, as we can 
see from Fig. Hfb), the TC modulus (solid line) goes to 
1 with a symmetric Lorentzian peak around the energy 
resonance, while the transmission phase (dashed line) 
smoothly changes by ir. In addition to the above res- 
onance two extremely sharp resonances for Ein around 
0.3 meV and 0.6 meV, are present. They are the so-called 
Fano resonances (with asymmetric line-shape), which 
can be ascribed to an effect brought about by electron- 
electron correlation 27 . In fact they originated from the 
interference between a resonant mechanism, due to the 
Coulomb blockade in the QD, and a non-resonant mech- 
anism, given by the transmission in the dot via a direct 
trajector y 23 ! 25 . In order to give a better insight into their 
properties we have reported in the top panel (c) a zoom 



of the modulus and phase of the TC at an initial energy 
of the incoming electron Ein around 0.3 meV. Here the 
modulus reaches 1 and then goes to zero, showing the 
typical asymmetric Fano lineshape. The transmission 
phase increases smoothly by 7r on the resonance peak, 
but shows an abrupt drop of 7r in correspondence of the 
zero of the transmission probability^ 7 .. We stress that 
the curves in Fig. [1] correspond to a single-channel two- 
particle scattering, with the energy of the transmitted 
electron fixed by the boundary conditions. As a conse- 
quence, no entanglement is generated between the two 
electrons here. The more interesting cases of two- and 
multi-channel scattering will be presented in the follow- 
ing section. 




FIG. 1: (Color online) Bottom panel (b): TC as a function 
of the initial energy of the incoming electron Ein- The solid 
line represents the modulus of the TC and refers to the left 
ordinate axis; the dashed line represents the phase of TC and 
refers to the right axis. Top panels: (a) profile of the potential 
V(x) in the scattering region: the potential well is 150 meV 
deep and 80 nm wide and is connected to the leads through 
two 10 nm barriers of 30 meV; (c) zoom of the TC modulus 
and phase for Ein — 0.3 meV where a Fano resonance occurs. 

In the last part of this section we describe the proce- 
dure we adopt to evaluate the entanglement. In fact as a 
consequence of the scattering, quantum correlations ap- 
pear between the energy levels E n of the electron bound 
in the potential well and the energies T n of the scattered 
electron allowed by energy conservation. Such an entan- 
glement may be evaluated from the transmitted compo- 
nent of the two-particle wavefunction in the right lead, 
given by Eq. ([3]). Obviously only the travelling compo- 
nents, with n < M in the right hand side of Eq. |3]), 
must be considered, since they are the only ones giving 
a non vanishing contribution to the current and also the 
only ones that could be revealed by measuring the en- 
ergy of the electron propagating in the right lead. On the 
other hand the evanescent component of the wavefunc- 
tion plays no direct role in the appearance of detectable 
quantum correlations. In this good entanglement 

measure is given by the von Neumann entropy of the re- 
duced density matrix pn e d, obtained by tracing the two- 
particle density matrix p = |** r )(* tr |/ {^ tr \^ tr } over 
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the degrees of freedom of one of the particles^, where 

1p(xi,X 2 ) tr = J2n=0 C nXn(x 2 )e lk "- Xl . 

pRed will be a (M+l) x (M+l) diagonal matrix defined 

as 



PRe 



diag 



I CO I 



cm 



with 



E 



M I I 
n=0 l C «l 



(4) 



(5) 



Thus the entanglement can be expressed by means of the 
von Neumann entropy as 



e = -Tr [pRed In pRed] = - 



M 

Ei 

n=0 



! ln 



(6) 



e is bound in the interval [0, In (M + 1)] with e = indi- 
cating no entanglement and s = In [M + 1) indicating a 
maximally- entangled state. We stress that, unlike other 
works estimating the quantum correlation in QDs 8 -^, we 
are not considering the transmitted and reflected compo- 
nent of the scattering wavefunction as two different states 
that can be entangled with the QD. We estimate the en- 
tanglement between the dot and the transmitted electron. 
In fact, in our approach the measure of the entanglement 
created in the system, evaluated by means of Eq. (|6|) , is 
not explicitely given in terms of the device transmission 
and reflection coefficients, but it is a function of the am- 
plitude probabilities of finding the transmitted particle 
in one of the possible energy states. Furthermore, it is 
worth noting that, although we do not use explicitely the 
criteria developed to treat the entanglement of identical 
particle o 5 ' 39 ! 40 , the electrons in our system are fully in- 
distinguishable. Nevertheless, since the energy of the in- 
coming electron is not sufficient to ionize the QD, we may 
safely assume that the scattered particle is far enough 
from the one left in a bound state of the dot so that the 
overlap between their spatial probabilities is negligible. 
In other words we use the spatial position of the elec- 
trons to "distinguish" them, while the entanglement is 
between their energy states, as suggested in Ref. HH. 



In panel (a) of Fig. [2] we report the entanglement e 
of the system as a function of the initial energy of the 
incoming electron and of the depth of the potential well, 
whose width is kept constant at L— 40 nm. The entan- 
glement presents extremely sharp peaks that, as we shall 
see, correspond to the resonances exhibited by the TCs 
of the various channels. From the 2D representation of 
the same data, displayed in Fig.^b), we observe that the 
maxima of the entanglement spread in the region corre- 
sponding to a three-channel scattering, that is separated 
by the dashed line from the one where the number of ac- 
tive channels is two. This suggests us that also the kind 
of process, i.e. two- or multi-channel scattering, plays an 
important role into the entanglement formation. The two 
cases will be better analyzed in the following subsections. 




FIG. 2: (Color online) Left panel (a): entanglement e as a 
function of the initial energy of the incoming electron Ein 
and of the depth of the potential well of 40 nm. Right panels: 
(b) 2D contour plot of the same data; the dashed line is a 
guide for the eyes and separates the zone corresponding to 
two-channel scattering (left) from the one corresponding to 
a three-channel scattering (right); (c) modulus of TC of the 
three channels as a function of Ein with h fixed at 100.5 meV. 
Note that the channel 2 is activated only at Ein — 21 meV. 



III. NUMERICAL RESULTS 

In this section we analyze the entanglement formation 
in a two- or multi-channel scattering, and its relation to 
the resonances in the TCs. In particular, it is of interest 
to study how the entanglement depends upon the kind 
of resonances, since the detection of the latter, usually 
feasible through a current-voltage characteristic of the 
quantum device, can give information about the former. 
In our approach we solve numerically the system for dif- 
ferent potentials V(x) (obtained by varying the depth 
h and the width L of the well), and for different initial 
energies of the incoming electron Ein- 



A. Two-Channel scattering 

Here we study the creation of the entanglement when 
the incident particle, as a consequence of the scattering 
with the particle bound in the QD, is transmitted with 
two possible energies, Tq and Xi, and, correspondingly, 
the hnal QD energy can be Eq and E\ . 

In the top panel of Fig. [3] we report the entanglement 
as a function of the kinetic energy of the incoming par- 
ticle, for the potential sketched in the inset of the figure. 
In particular we consider a width of the potential well 
of 30 nm and a depth of 110 meV. At low energies the 
scattering does not lead to the appearance of quantum 
correlations between the two particles, since only a single 
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channel is possible for the transmission. It is therefore 
possible to attribute a specific energy to each particle: 
To for the scattered electron and Eq for the bound one. 
When Ejn reaches a threshold value of about 14 meV a 
new channel comes into play as it can be seen from the 
bottom panel of Fig. [3j There, the dependence of the 
modulus of the TC is reported as a function of the initial 
energy of the incoming electron. In correspondence to 
the energy at which the second channel is activated, the 
entanglement shows a sharp increase. 
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FIG. 3: (Color online) Top panel: the entanglement e as a 
function of the initial energy of the incoming electron Ein- 
The inset displays the potential V(x) inside the scattering 
region: we consider a potential well 110 meV deep and 30 nm 
wide. Bottom panel: the modulus of the TC of the channels as 
a function of the initial kinetic energy of the incident electron: 
0-th channel (solid line) and 1-st channel (dashed line). The 
scale of the abscissa is the same in the two panels. 



By comparing the two panels in Fig. [3j it is clear that 
the resonances of the TCs play a key role in the entan- 
glement formation. In particular, the behavior shown by 
the entanglement in correspondence of a Fano resonance 
is very different from the one exhibited at a Breit-Wigner 
resonance. When Ein is around 15 meV, a Fano reso- 
nance of the TC is observed for channel 0. To better 
understand this phenomenon, we report in the top panel 
of Fig. [3] a zoom of the curves of modulus and phase (the 
latter was not shown in Fig. [3J for clarity) of the TCs of 
the two channels. It is worth noting that, unlike the case 
of the single-channel scattering, here the modulus of the 
TC attains small but non zero values before the max- 
imum, which, in turn, results to be significantly lower 
than 1. The TC of channel 1 shows, in correspondence 
of the Fano resonance of channel 0, a Breit-Wigner res- 
onance characterized by a phase change of about 7r. In 
the energy interval around 15 meV the behavior of the 
quantum correlations, appearing in the systems as a con- 
sequence of the scattering, is peculiar (see the top panel 
of Fig. 0]). In fact when the modulus of the TC of the 
channel reaches its lowest value, the entanglement curve 
presents a minimum. Such a minimum is placed between 



two very close maxima, where the entanglement (evalu- 
ated, as usual, by means of the von Neumann entropy 
of the reduced density matrix) is equal to In 2. This 
value indicates the condition of maximal entanglement 
in a two-channel scattering. Such a condition is reached 
when | Co | 2 and |ci| 2 , i.e. the probabilities that the scat- 
tering occurs through the channel or 1, respectively, are 
equal, and it implies that the lack of knowledge about the 
state in the one-particle subspace is maximum. We also 
report, in the top panel of the Fig. [5] (dashed line), the 
ratio of the two transmission probabilities: the entangle- 
ment is maximum when the scattering probabilities in 
the two channels are the same, as indicated by the hor- 
izontal dotted line drawn as a guide for the eyes. What 
we found here is in agreement with previous analyses on 
the two-electron entanglement production in two-electron 
systems for a two-channel scattering model, where both 
particles are injected in only one of two leads&2£. In fact, 
also in those cases, the entanglement shows a maximum 
when the transmission probabilities for the two channels 
are identical, while it vanishes in correspondence of the 
single-particle resonances, where there is no uncertainty 
about the energy of the particle. 
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FIG. 4: (Color online) Top panel: the entanglement e (solid 
line) and ratio between the probabilities of finding the scat- 
tered particle in the transmission channel 1, |ci| 2 , and 0, |co| 2 
(dashed line), against the initial energy of the incoming elec- 
tron Ein, close to a resonant condition. Note that, while e 
can range from (no entanglement) to In 2 (maximally en- 
tangled state), the bottom of the left scale is set to 0.25 to 
optimize the magnification effect in the figure. The horizontal 
dotted line is a guide for the eyes indicating the point of equal 
scattering probabilities for the two channels, i.e. 1 in the left 
ordinates axis. Bottom panel: TC of the two channels as a 
function of Ein around a Fano resonance. The solid and the 
dash-dotted lines represent the modulus of TC of the channel 
and 1 respectively, and refer to the left ordinate axis; the 
dashed and the dotted lines represent the phase of the TC of 
channel and 1, respectively, and refer to the right axis. 



Figure [3J shows that the TC of the channel presents 
a Breit-Wigner resonance for Ein — 17.5 meV. Even if 
the modulus of the TC becomes equal to 1, the entan- 
glement curve does not display maxima or minima. This 
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is due to the fact that the Breit-Wigner resonance of the 
TC of channel does not influence the TC of channel 1 , 
which does not show, at the specific energy, resonances 
of any kind. Therefore, the scattering phenomena taking 
place in the energy interval around 17.5 meV do not play 
a special role into entanglement formation of the two- 
particle system. Such a behavior is in agreement with 
the one observed in other works, showing that the maxi- 
mal value of the conductance does not always correspond 
to the maximal entanglement 3 ^. 



B. Multi-channel scattering 

Let us consider now, the case of a multi-channel scat- 
tering (the scattered particle can leave the QD with more 
than two energies), and let us investigate the different 
roles played by the Fano and Breit-Wigner resonances in 
the entanglement formation. 

Here, we consider a potential well 150 meV deep and 
40 nm wide, as reported in the inset of Fig. \5\ For Ein 
around 21 meV, the scattering passes from a two-channel 
to a three-channel process, and this transition is charac- 
terized by a sharp increase of the entanglement. The 
same behavior is found at Ein = 40 meV, where an ad- 
ditional transmission channel is switched on (see Fig. [5]). 
This is in agreement with the results of the previous sec- 
tion and may be considered representative of a more gen- 
eral behavior, occurring whenever a new channel becomes 
effective in the scattering process. 




FIG. 5: (Color online) Top panel: the entanglement e as a 
function of the initial kinetic energy of the incoming electron 
Ein- The inset shows the profile of the potential V(x) in 
the scattering region: the potential well is 150 meV deep and 
40 nm wide. Bottom panel: modulus of the TC of the four 
channels as a function of Ein'- channel (solid line), 1 (dotted 
line), 2 (dashed line) and 3 (dash-dotted line). 



The quantum correlations appearing in the system, as 
a consequence of a multi-channel scattering, for the en- 
ergy values around a Fano resonance show some differ- 
ences with respect to the two-channel process. In Fig. [B] 



we report the curves of the entanglement and, in the 
insets, the modulus and phase of the TCs of the three 
channels (0, 1, 2 from top to bottom) against Ein, in 
the energy region around Ein — 27.5 meV, where a Fano 
resonance occurs for channel 0. In fact, we observe, from 
the uppermost inset that also in this case the modulus of 
the TC has a local minimum before reaching the maxi- 
mum. The type of resonances of channels 1 and 2 (mid- 
dle and bottom insets) cannot be clearly identified by the 
module of their transmission coefficients, whose peaks are 
almost symmetric. However, it is clear for the transmis- 
sion phases, that channel 1 exhibits a Fano resonance, 
with the phase essentially unchanged through the peak, 
while channel 2 shows a Breit-Wigner resonance, with a 
global phase variation of w. We stress that, unlike the 
two-channel scattering case, here the entanglement does 
not show a minimum. Such a behavior can be ascribed to 
the fact that, when the modulus of the TC of channel 
is small, the TCs of the other two channels attain values 
comparable to each other. This means that the probabili- 
ties of finding the scattered particles in those channels are 
almost equal, and there is still a lack of knowledge about 
the state of the one-particle subsystem. Furthermore, we 
note that here the entanglement presents a single maxi- 
mum whose value exceeds In 2. Actually, the fact that the 
number of degrees of freedom is larger than 2 increases 
the uncertainty about the constituents of the system; in 
fact Eq. ([6]) gives a maximum value for the amount of 
quantum correlations that is larger for a larger number 
of possible states. For example, a system of two qutrits 
is able to attain a larger value of e than a two-qubit 
system, — The ability to tailor not only the degree of the 
entanglement, but also the number of possible states of 
the two subsystems, by independently tuning Ein and 
the QD confining potential, could also have implications 
beyond the theoretical estimation of the entanglement. 
The behavior described above is repeated at Ein — 50 
meV where two Fano resonances occur for the second and 
third channel in a four-channel scattering (Fig. 0. ft is 
worth noting that, in correspondence of a Breit-Wigner 
resonance in the channel for Ein = 25 meV, no ad- 
ditional resonance occur in the other channels and the 
entanglement does not present maxima or minima as it 
can be clearly seen from the upper panel of Fig. There- 
fore also in the case of multi-channel processes, the Breit 
Wigner resonance seems not to induce sharp variations 
of e. 



IV. CONCLUSIONS 

The controlled production and detection of entangled 
particles in the solid state environment represents an 
experimental challenge. In this spirit, various propos- 
als for producing bipartite entangled fermionic systems 
have been advanced, on the basis of different physi- 
cal mechanisms requiring a direct interaction between 
particles^£&£ ' 10 i 33 . In this paper we have investigated 
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FIG. 6: (Color online) The entanglement e as a function 
of Ein around a Fano resonance of the TC of channel Oin 
a three-channel scattering. The insets display the modulus 
(solid line, left axis) and phase (dashed line, right axis) of TC 
of three channels (0, 1, 2 from top to bottom) in the same 
energy range of the main graph. 



the quantum correlations appearing, as a consequence of 
a Coulomb scattering, between two electrons having the 
same spin, in a system of physical interest, where the de- 
gree of the entanglement results to be controllable by a 
proper tuning of the carrier energy and of the QD poten- 
tial. Such a system consists of a quasi ID double-barrier 
resonant tunnelling device, where an electron incoming 
from one lead is scattered by the potential structure and, 
via the Coulomb interaction, by another electron bound 
in the QD. 

The numerical procedure, used to solve the model, is 
a generalization of the quantum transmitting boundary 
metho d 25 i 36 . It permits us to obtain the reflection and 
transmission amplitudes of each scattering channel, for 
various configurations of the potential, as a function of 
the initial energy of the incoming electron. However, we 
stress that, unlike the approaches followed by Lopez et 
alM and by Oliver et here we did not use the re- 
flected component of the scattered electron wavefunction 
to evaluate the entanglement of the two-particle system, 
but we estimated the quantum correlations showing up 
between the QD eigenstates and the transmitted parts 
of the electron wavefunction. Furthermore, this proce- 
dure makes it possible to investigate the role played by 
the resonances of the transmission spectra into the en- 
tanglement. Although our numerical analysis has been 
performed by using the GaAs material parameters, they 
can be considered representative of a more general be- 
havior. 

Our simulations show that the entanglement depends 
upon the kind of the resonance appearing in the transmis- 



sion spectrum. A single Breit-Wigner resonance is found 
not to induce peculiar effects on quantum correlations. 
On the contrary, in correspondence of a Fano resonance 
of one of the TCs, not only the other channels TCs ex- 
hibit local maxima, but also the entanglement presents 
sharp peaks. Such a behavior can be related to the na- 
ture of the resonances themselves. In fact a Breit-Wigner 
resonance is essentially a one-particle effect, showing up 
also in a single-electron scattering. On the other hand, a 
Fano resonance is a genuine multi-particle phenomenon, 
due to electron-electron correlation, and is not present 
for single-particle systems. 

Furthermore, we showed that the appearance of quan- 
tum correlations in our system is also affected by the 
number of the transmission channels, i.e. the number 
of possible energy levels of the scattered particle (and, 
due to the energy conservation, of the bound particle). 
In fact, the entanglement shows a sharp increase when- 
ever a new channel is turned on. Moreover, its behavior 
for energy values around a Fano resonance is found to 
depend upon the kind of process: two- or multi-channel 
scattering. For the two-channel case, the entanglement 
presents a minimum between two close maxima, which 
indicate the maximal uncertainty about the state of the 
system. In the multi-channel case, a single maximum of 
the entanglement, with no minima, is observed. When 
the energy levels of the scattered and bound electrons 
are only two, the minimum of the entanglement is found 
in correspondence of the local minimum of the TC of 
the Fano resonant channel. In this case it maximizes the 
possibility to ascribe specific energy states to the subsys- 
tems. On the other hand, for a multi-channel scattering, 
a minimum of TC of a Fano resonant channel does not 
imply a decrease of uncertainty about the subsystems, 
since the TCs of the other non Fano resonant channels 
attain values comparable to each other. 

Finally, the results of our paper suggest that the ma- 
nipulation of Fano resonances and of the number of 
scattering channels may allow to significatively influ- 
ence the degree of entanglement between the transmitted 
electrons and the QD. A promising development of the 
present work could be the study of the entanglement in 
the case of scattering of a single electron by a few charges 
confined in the QD, in connection with experimental re- 
sults obtained for the coherent components of the trans- 
mitted current in the case of the multi-occupancy of the 
do t 18 ' 42 . In the latter case, with three or more particles in 
the system, the spin degrees of freedom cannot be factor- 
ized and their inclusion in our approach, although quite 
feasible, results very challenging form the computational 
point of view. 
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